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ABSTRACT 

We performed a detailed analysis of extensive photometric observations of a sample of 
most active dwarf novae, that is SU UMa stars which are characterised by supercycle 
lengths shorter than 120 days. We found the observational evidence that supercycle 
lengths for these objects have been constantly increasing over the past decades, which 
indicates that their mean mass transfer rates have been decreasing during that time. 
This seems to be a common feature for this type of stars. We present numerical results 
in each case and estimate time scales of future development of these systems. This 
study is important in the context of evolution of dwarf novae stars and perhaps other 
cataclysmic variables. 
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1 INTRODUCTION 



As it was pointed out in one of the papers of IPaczvriskH 
i|l997l) . massive professional photometric surveys are to- 
day sources of m any, often unexpe cted, discoveries (e.g. the 
ASAS project of Poimanskl fl997) or the Kepler mission of 
iBorucki et all (|2010D ). On the other hand, databases collect- 
ing amateur photometric measurements are growing rapidly, 
as tel escopes and CCD detectors are easily accessible nowa - 
days (|Pattersonl l201ll : ISzkodv et all 120121 : IWaagenl |2012| ). 
This gives us access to a number of valuable open astronom- 
ical archives, both professional and amateur, with a richness 
of data gathered during the past decades. We have now a 
great opportunity to trace photometric behaviour of a lot of 
variable objects, for instance dwarf novae stars. 

Dwarf novae are cataclysmic binaries with a white dwarf 
as a primary and a low mass main sequence star, filling its 
Roche lobe, as a secondary. Due to the transfer of mass 
from the secondary onto the non-magnetic primary, an ac- 
cretion disk is formed. Such a mass transfer leads to out- 
bursts in dwarf novae, as they are caused by sudden gravi- 
tational energy release due to accretion of this material onto 
the white dwarf. Dwarf novae with shortest orbital periods 
(Port < 2.5 h) are called SU UMa stars. They are char- 
acterised by the presence of superoutbursts, next to nor- 
mal outbursts, which are about 1 mag bright er and about 
ten times less frequent than normal outbursts (|Hellierll200il : 
I Warner! [20031 ). 

The length of the supercycle, i.e. the time between two 
successive superoutbursts, is one of the most fundamental 
properties of SU UMa stars and is specific for each of them. 
The nature of this feature is still uncertain. Nowadays there 
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is an open debate about the cause of superoutbu rsts. The 
thermal-tidal instability (TTI) model (|Osakill 19961 ) suggests 
that the sypercycle length (P S e) is set by the mass transfer 
rate (Mtr)- According to this scenario, the observed P sc is 
inversely proportional to Mtr- The TTI model shows that 
normal outbursts remove less material from the disk than 
it has been accreted since the last normal outbursts. This 
makes the disk grow during a number of normal outbursts. 
When the disk reaches the critical radius (r ~ 0.46a, with a 
being the binary separation), it becomes eccentric and tidal 
effects occur, which causes superoutburst, during which the 
disk is shrinking again. Although the TTI model is com- 
monly accepted, it still seems to need some improvements. 
For instance the TTI simulations which reproduced the light 
curve of one of the most active dwarf novae, ER UMa star, 
required an artificial increase of the mass transfer rate by 
the factor of ten in comparison to values expected fro m the 
theory based on gravitational radiation (Osaki 1995). An- 
other example is the fact shown by ICannizz 

3 (|2012L and 

references therein) that embedded precursors of superout- 
bursts are present also in systems with long orbital peri- 
ods. This argues for a more general model, which is not 
restricte d to the mass ratio o f q < 0.25, as the TTI model. 
Recentlv lOsaki fc Katol l|2012l ) presented evidence in favor of 
this model. Although their analysis was based on variations 
of the negative superhump period of a single dwarf nova, 
V1504 Cyg, they claimed th at the cause of superoutbursts 
is finally revealed. However, ISmakl (|20 13f ) showed that this 
object cannot be considered as a representative for all sys- 
tems which show negative superhumps. He also presented 
a number of arguments which suggest that superoubursts 
are caused b y strongly enhanced mass transfer (EMT) rate 

i|Smakll200Sl ). 

Recent observations l|Otulakowska-Hypka et all 120131 ) 
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surprisingly showed that the supercycle length of one of the 
most active dwarf novae, IX Draconis, has been increasing 
with a constant rate since the last twenty years. This is 
interesting in the context of evolution of such systems. En- 
couraged by this fact, we decided to investigate the same 
issue for other well-observed active objects of the SU UMa 
type. 

The paper is arranged as follows. In Section [2] we give 
information on the data used in this research. The analysis 
is presented in Section [3] Results together with discussion 
are given in Section [4] We summarize the main conclusions 
in Section [5] 



2 DATA 

We searched accessible archives of amateur astronomical ob- 
servers and automatic professional surveys to create joint 
and as complete light curves as possible. In a few cases we 
also have an additional data from our previous observational 
campaigns and from other observers. The relevant photo- 
metric data was obtained from th e following databases: the 
project l|Poimahskil [l997j h the MEDUZA0 project 
of the Variable Star and Exoplanet Sect ion of the Czech 
Astro nomical Society (|Brat fc Zeidall2010t) . the Kepler mis- 
SIO rFI (|Borucki et al.ll2010l ). and the databases with amateur 
observations: AAVS CQ, AFOEV0, and BAAVSS0. 

At first, we intended to perform the analysis for most 
of the well-observed (which often means the brightest) ac- 
tive SU UMa stars. However, we found not enough data for 
the following objects: BF Ara, BK Lyn, BR Lup, CI UMa, 
IX Dra, MN Dra, NY Ser, SDSS J210014, and UV Gem, to 
reach any unambiguous conclusions. We selected only those 
targets which have really good time coverage of their light 
curves, i.e. those, for which we were able to easily distinguish 
superoutbursts from the rest of the signal. The chosen ob- 
jects have light curves spreading over at least twenty years. 
In Table [1] we present the final list of objects together with 
information on the used databases for each star. 

Although our selected objects are all of the SU UMa 
type and they are located below the period gap in the or- 
bital period distribution for dwarf novae, they are diverse in 
terms of the evolutionary status. Figure Q] shows the relation 
between the orbital period-superhump period excess (e) and 
th e orbital period (P o r b) for dwarf novae stars, noticed first 
bv lStolz fc Schoembsl i|l984h . where we marked by squares 
positions of the selected objects. Thanks to the dependence 
between e and the mass ra tio (q) for dwarf novae stars re- 
ported bv lPattersonl (119981 ): 



0.23g 
l + 0.27q' 



(1) 



we can use the relation between P or b and e as an excellent 
plane to examine evolution of these stars, since the mass 
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Figure 1. Relation between the orbital period-superhump period 
excess (e) and the orbital period (P or b) for dwarf novae stars. 
The objects used in this research are shown by squares. The line 
indicates possible positions of RZ LMi star. 



ratio decreases with time due to the mass loss from the sec- 
ond ary. In Fig.pjsmall p oints are known dwarf novae objects, 
from lOlech et alj l|2011f ). It was impossible to mark the exact 
position of one of the stars, RZ LMi, since its orbital period 
is unknown. However, based on the value of its superhump 
period, we know that it must be located somewhere on the 
solid line, which is shown in this figure. 



3 ANALYSIS 

To create final joint light curves from all data sources listed 
in Tab. [1] we decided to use only the certain data points, 
without upper limits or visual observations. For all the data 
sources, with the exception of the Kepler data, the merger of 
light curves was simple and straightforward. The only thing 
to do was to bring all the measurements into the common 
final time system, the HJD. For light curves from the Kepler 
mission this process was different. Here we had to transform 
the flux into magnitudes. For this we used the Pogson's equa- 
tion plus an arbitrary chosen constant shift, since we had no 
reference point. We did this step very carefully, but it is 
still possible that there exists some error between magni- 
tudes from the Kepler part of the data and the rest of data 
sources. Even if this was true, it would be insubstantial, since 
our study concerns time series analysis, not the brightness. 
In the case of the Kepler data the time was originally given 
as BJD-2454833.0. The rest of the data was in the HJD. We 
did not have to convert one time system to another because 
the difference between these two sta ndards is of the order of 
a few seconds l|Eastman et alj|2010h . which is a few orders 
of magnitude below our precision of determination P ac . 

As a result we obtained a huge data sets with a number 
of photometric measurements of the order of tens of thou- 
sands for each of the objects. 

In Fig. [2] we show as an example the final light curve 
for YZ Cnc star, which was used in the further analysis. 
The obtained light curves of all the objects are not uniform. 
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Table 1. List of objects with information on the used databases in this study. 



Object 


ASAS 


MEDUZA 


AAVSO 


AFOEV 


BAAVSS 


Kepler 


extra data 


DI UMa 




X 


X 


X 


X 




Rutkowski et al. (2009) 


ER UMa 




X 


X 


X 


X 




Kato & Kuniava (1995) 


RZ LMi 




X 


X 


X 


X 




Olech et al. (2008) 


SS UMi 




X 


X 


X 


X 




Olech et al. (2006) 


V1159 Ori 


X 


X 


X 


X 


X 






V1504 Cyg 




X 


X 


X 


X 


X 




V344 Lyr 






X 


X 


X 


X 




V503 Cyg 




X 


X 


X 


X 






YZ Cnc 


X 


X 


X 


X 


X 
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Figure 2. Example of a final light curve of YZ Cnc used in the further analysis (top). Zoom in the most recent part of this light curve 
where superoutbursts are clearly visible (bottom). 



There is a yearly trend present in all of them, with the excep- 
tion of the Kepler data, caused by the variable observational 
conditions. In addition, the density of observations increases 
with time, as the number of amateur observers is constantly 
growing. 

For this unev enly sampled data, we decide d to use the 
ANOVA software |Schwarzenberg-Czerny|[l996l ) for the time 
series analysis. First, we divided each of these light curves 
into five bins of equal time range. For each of the bins we 
searched for the most prominent peak in the corresponding 
power spectrum obtained from ANOVA. The data suffered 
from aliasing because of the under-sampling. However, we 
were able to derive values of P sc for vast majority of the 
time bins for all of the objects. These values correspond to 
the most dominant peaks in a reasonable range of frequen- 
cies of periodograms. An example of such power spectrum 
is presented in Fig. [3] 

The obtained value of supercycle length for a given time 
bin is naturally not exactly equal to each time period be- 
tween each pair of successive superoutbursts during the time 
covered by the time bin, because not all of them are strictly 
regular. There are some subtle differences between the P sc 
measured on short time scales, as w as already shown many 
times in th e literature, for insta nce bv lAntonvuk fc Pavlenkcl 
(2005) and lZemko et~al1 (|2012T ). However, our analysis aimed 
at examination of the overall behaviour of the P 3C , i.e. on 



40000 

| 30000 

™ 20000 
> 

S 10000 
< 













































-4 
























\ _ 











0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 
Frequency [c/d] 



Figure 3. An example of a resulting periodogram. This is the 
case for the last time bin of ER UMa star. The most prominent 
peak here corresponds to the P sc = 50.9 ± 0.5 d. 

long time scales of the order of decades. Thus, based on 
the ANOVA statistics, we derived one value of P sc for each 
time bin, which corresponds to the range of observation time 
from a few years to one decade, depending on the object. All 
the results are presented in the next section, and the exact 
ranges of time bins for all stars are given in Table [2] 



4 RESULTS AND DISCUSSION 

In Fig.[4]and[5]we show the obtained evolution of supercycle 
lengths for all of the analysed objects. Unfortunately, the 
quality of the rest of the data did not allow us to investigate 
more objects, so the statistics is poor. Seven out of nine 
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Figure 4. Results from ANOVA which show increasing supercycle lengths. Data points without uncertainties are possible but not certain 
due to poor light curves coverage. 
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Figure 5. Results from ANOVA which show constant supercycle lengths. Due to the light curves quality, for some of the time ranges 
the analysis was impossible. We show only the available measurements. 



examined targets show clearly growing supercycle lengths 
during the past decades (Fig. \$$. The only two examples 
which seems to have a constant value of P sc are the ones with 
the poorest time coverage of the observations, before the 
Kepler data, which could have influenced the results (Fig.[SJ). 

The behaviour of the P sc of YZ Cnc (the last one in 
Fig. [3} is a puzzle. First of all, it has the highest value of its 
period among all of the objects. Secondly, there is a peculiar 
drop in the P sc between the third and the fourth time bin, 



which looks like a transition from the increasing to constant 
supercycle length. It is possible that this is only a short time 
scale fluctuation, like for instance in the case of V503 Cyg. 
However, this drop is significant in comparison to the un- 
certainties of a single measurement, unlike in the case of all 
the rest of the objects, which have fluctuations of the order 
which is about one magnitude smaller. Since this effect is 
so strange and unique, we checked it twice, and we claim 
that this is certainly not an artificial result. Because of the 
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Figure 6. Plot from Otulakowska-Hyp ka et al. (2013) which 
shows the increasing supercycle length of IX Dra during the past 
twenty years. The data points are taken from the literature. Un- 
certainties are given when available. The line represents the best 
fit to the data. The corresponding rate of increase of the period 
is P = 1.8 X 1CT 3 . 



Table 2. Rates of changes of P BC which correspond to the best 
fit to our measurements, presented in Fig. [3] [3] and \E\ together 
with respective ranges of P sc and HJD. 



Object Psc P™ in P™ ax HJD range 

[days] [days] [years] 



fact that the value of the supercycle length for YZ Cnc is 
the highest one, we cannot exclude a possibility that there 
is some critical value for such a high P sc which could have 
caused the drop. 

For completeness, in Fig. [6] we additionally 
present the plot from ou r previous publication 
l|Otulakowska-Hvpka et all I2013T ). which shows the in- 
creasing supercycle length of another active ER UMa-type 
star, IX Draconis. This plot indicates that the rate of the 
period change is P — 1.8 x 10~ 3 , see FigdJ This result 
is based on four independ ent measurements f rom various 
papers from the l i terature ([I shioka et al. 2001; lKloselll995l ; 
lOlech et all 12004 lOtulakowska-Hypka et all l2013h . which 
were obtained from short observational campaigns in each 
case. This star was in our original list of objects for the 
present study, however, the examination of its whole light 
curve was impossible since it was too poor to perform the 
ANOVA analysis. 

For all of our targets we found positive values of their 
period derivatives of supercycle lengths, P sc . They are pre- 
sented in Table [2] 

Only in two cases, i.e. V344 Lyr and V1504 Cyg, the P 3C 



values are so small that they seems to be constant. Within 
the uncertainties they could also be either positive or nega- 
tive, thus the conclusive distinction is impossible here. For 
two other objects, DI UMa and V503 Cyg, it is also possible 
that within the uncertainties the values of P sc are positive, 
constatnt, or negative. This is caused by the fact that mea- 
surements of their P sc show the highest dispersion. Nonethe- 
less, for these two stars growing trends of their supercycle 
lengths are clearly visible in the Fig. U 

A comparison of these results with other examples 
known from the literature is satisfactory. They are in 
agreement at least with the order of magnitude. In the 
case of the ER UMa star, we found two previously pub- 
lished val ues of its supercycle length derivative: P 3C ~ 
4 x 10~ 3 (iRobertson et al.lll995h . and P sc = 6.7(6) x 10~ 4 
(|Zemko et al.|[2012l ). The first one is based on observations 
which were gathered between 1993-1995 (HJD « 2449000), 
so it is incomparable with our result. However, authors of 
this publication show that the mean supercycle length equals 
42.95 ± 0.06 d, which is in a perfect agreement with the P sc 
for HJD w 2449000 in our Fig. H On the oth er hand, the 
second result presented by IZemko et al.l l|2012h is based on 
analysis of a light curve with the time coverage of twenty 
years. Here the result is consistent with ours to the order 
of magnitude, and the small difference probably comes from 
the differences in the data sources and methods of analysis. 
Another obje ct, RZ LMi, has also a p ublished value of Psc = 
-1.7 x 10" 3 l|Robertson et al.lll995T ). This one is based on 
observations taken between 1992-1995 and is completely dis- 
crepant with our result. This is probably caused by the fact 
that this observing time was not long enough to reach any 
other conclusion, and it perha ps reflects some mi nor fluctu- 
ations of the supercycle length, fkato et al. (2000) found the 
mean supercycle for SS UMi to be equal to 84.7 d for the 
observations time HJD ~ 2451000, which is again in a per- 
fect agreement with our Fig. [3J The same is for V1159 Ori, 
with P 3C = 45 d ar ound HJ D w 2450000, as presented by 
ISzkodv et all (|l999T ). In turn lAntonvuk fc Pavlenkol (|2005T l 
showed for the five years spread of observations that the su- 
percycle for V1504 Cyg seems to stay more or less constant, 
apart from short time scale fluctuations. 

Our examination of supercycle lengths' behaviour on 
long time scales demonstrates that all of the analysed active 
SU UMas have a positive value of their superoutbursts pe- 
riod derivative, and eight out of ten objects have a very 
clearly increasing trend in their P sc evolution. This be- 
haviour seems to be general for this type of stars. The fact 
that for most of these sources there is a clear evidence of 
increasing supercycle lengths means that their mean mass 
transfer rates has been constantly decreasing during the 
past decades. This scenario of evolution is in agreement 
with results for extreme ly interesting ER UMa -type object, 
BK Lyncis, presented bv lPatterson et alj l|2012t ). In this pub- 
lication the authors showed a hypothesis of the evolutionary 
path of this star as its ER UMa stage is a transient phase of 
evolution, preceded by classical nova and novalike phases. If 
we assumed that this is true for all active SU UMas of our 
sample, this would mean that these objects are constantly 
fading since their classical nova eruptions. Considering the 
fact that they are rather diverse in terms of their orbital pe- 
riods, i.e. evolutionary stages (see Fig.[T}, even though they 
are all below the period gap, we can try to predict time 
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29 


V1159 Ori 
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-4 


42.6 


55.1 


19 


V1504 Cyg 


(0.3 ±5.4) x 10" 


4 


110.6 


120.6 


20 


V344 Lyr 


(0.2 ± 2.5) X 1Q- 


4 


112.9 


123.9 


30 


V503 Cyg 


(0.6 ± 1.7) x 10" 


3 


112.1 


128.1 


20 


YZ Cnc 


(30.2 ± 1.4) x 10 


-4 


90.5 


149.5 
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Table 3. Estimates of time scales of next steps of evolution of 
our objects under the assumption of the constant growth of their 
supercycle lengths. The second column shows in how many years 
a star will become a typical SU UMa star with P ac = 300 d. 
The third column shows in how many years period bouncers will 
attain the stage of WZ Sge star with P ac = 10 yr. 



Object SU UMa stage WZ Sge stage 
(P ac = 300 d) (P ac = 10 yr) 



DI UMa 


1700 


23000 


ER UMa 


500 




IX Dra 


400 


6000 


RZ LMi 


1500 




SS UMi 


200 




V1159 Ori 


600 




V1504 Cyg 


20000 




V344 Lyr 


33000 




V503 Cyg 


800 




YZ Cnc 


100 
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scales of their further steps of evolution. All of our objects 
are extremely active as for the SU UMa class. The most ac- 
tive of them are sometimes distinguished as a separate class 
of ER UMa-type, with P sc m 50 d. They are believed to 
be typical SU UMas with only higher mass transfer rates 
(DI UMa, ER UMa, IX Dra, RZ LMi, and V1159 Ori from 
our sample). The rest of the sample are representatives of 
the SU UMa class of stars but still with extremely short su- 
percycles (SS UMi, V1504 Cyg, V344 Lyr, V503 Cyg, and 
YZ Cnc) . A typical supercycle length for a star of SU UMa 
type is equal to a few hundr ed days and f or the WZ Sge class 
it is of the order of decades (|Hellieril200ll ) . We adopted these 
standard values to find the expected time scales in which our 
objects will reach next levels of evolution. For all of the stars 
we estimate the time which is needed to become SU UMa- 
type object with P ac = 300 d. For the two period bouncers in 
our sample, DI UMa and IX Dra, we also give times needed 
to reach the next step that is becoming a WZ Sge star with 
P ac = 10 yr. Results of this estimate are shown in Table [3] 



5 CONCLUSIONS 

We analysed a set of photometric data covering last decades 
of observations for the most active SU UMa stars in order to 
study changes of their supercycle lengths. Despite the com- 
mon feature of the selected objects, which is the highest ac- 
tivity in this class of stars, they are diverse in terms of their 
orbital periods below the period gap. For all of the anal- 
ysed objects we found positive values of period derivatives 
of their supercycle lengths. There are some subtle fluctua- 
tions for short time scales in the behaviour of P ac , but the 
general trend is the same in each case. Increasing supercycle 
lengths mean that the mean mass transfer rates have been 
decreasing for these objects over the last decades. This is 
in agreement with the scenario of th e evolution of BK Lyn 
presented by iPatterson et all (|2012l ). which seems to be a 
general case. This phenomenon is important in the context 
of evolution of such systems. 
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